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NOMENCLATURE 


Symbol 


d 

G 

LWC 


’'i 

V 

w. 


Greek 

Symbol 

a 

6 


DeClnltlon 

Volume median water droplet diameter 

Glycol mass fraction in glycol-water solution 

Liquid water content in atmosphere 

Mass rate of water droplet impact per unit area on wing 
surface 

Ambient freestream temperature 

Freezing temperature of glycol-water solution 

Stagnation temperature 

True airspeed 

Equivalent airspeed 

Rate of flow of glycol solution per unit area 

Glycol mass fraction in glycol-water solution pumped 
through porous leading edge 


Angle of attack 
Local catch efficiency 


Subscripts 

L Lower panel 
U Upper panel 
max Maximum 
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INTRODUCTION 


The concept of pumping a glycol-water solution through a porous 
leading edge skin to achieve ice protection is not new. Operational 
systems using this concept have been employed on European airplanes 
for many years. However, no U.S. manufacturers have used liquid anti- 
ice protection for wing and tail leading edges up to this time. Re- 
cently there has been a surge of interest in the U.S. in this concept 
because of the potential advantages it offers in comparison with pneu- 
matic boots and hot air bleed systems. 

Several icing tunnel tests with liquid ice protection systems have 
been conducted in the NASA Lewis Icing Research Tunnel during the past 
few years to add to a relatively meager data base. The purpose of this 
report is to present a proposed method of analytically predicting the 
minimum fluid flow rate required to provide anti-ice protection with a 
porous leading edge system on a wing under a given set of flight condi- 
tions. Results of the proposed method are compared with the actual re- 
sults of an icing test of a real wing section in the NASA Lewis Icing 
Research Tunnel. This work was conducted under NASA Grant NAG 3 - 71 . 

SYSTEM DESCRIPTION 

A fluid ice protection system consists of a porous leading edge 
panel attached to the leading edge of a wing and a pump that dis- 
tributes a glycol-water solution from a reservoir to the leading edge 
panel through nylon tubing. The fluid flows through the panel onto the 
surface of the wing, providing either an anti-icing capability by dis- 
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solving supercooled water droplets and preventing the formation of ice, 
or a de-icing capability by chemically breaking the bond of established 
ice. A significant feature of the system is that protection is obtained 
aft of the panel by the flow of the fluid along the chord to the trail- 
ing edge, which often prevents the formation of ice anywhere aft of the 
active leading edge. 

The porous panel material most commonly used consists of two or 
three layers of stainless steel wire cloth that are rolled, sintered, 
and finish rolled to proper thickness . Recent development programs 
have also produced porous panels made of laser-drilled titanium sheet 
and various composite materials. 

A typical cross section of a porous leading edge panel installa- 
tion is shown in Figure 1. The edges of the active portion of the 
panel must be placed such that extreme positions of the stagnation 
point for which icing protection is required are not too close to the 
edge to prevent fluid from being distributed on both the upper and 
lower surfaces of the wing. 


PREDICTION METHOD 


Anti-ice protection is obtained by providing a glycol-water solution 

on the leading edge of the wing that mixes with the atmospheric water 
droplets as they impact on the leading edge. To provide anti-ice pro- 
tection, the resulting solution must have a glycol mass fraction that 
is high enough to prevent any freezing at the leading edge or on the 
wing upper and lower surfaces as the fluid flows aft along the wing. 
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Figure 1: Cross Section of a Porous Panel Installed 

on a Wing Leading Edge. 




The water droplets do not impact as a uniform mass rate on the lead- 
ing edge. The mass rate per unit area tends to be highest at or near the 
stagnation point and decreases in a chordwise direction on either side of 
the stagnation point. The distribution of water mass rate on the leading 
edge is described by the nondiraensional catch efficiency parameter 3. 

The acttjal local water mass rate of impact per unit area is given as 

M = 3(LWC)(V). (1) 

w 

A typical distribution of 3 is shown in Figure 2. The trajectories of 
the water droplets relative to the flow streamlines determine the 3 
distribution. Thus 3 is a function of the airfoil shape and size, the 
remote airspeed, the air density, and water droplet diameter. Given 3, 
and the airspeed, the local and total mass rate will be directly propor- 
tional to liquid water content. 

Until recently, the only methods of predicting 3 were empirical 
and semi-empirical. However, the development of improved methods of 
computational aerodynamics has resulted in several different computer 
codes that predict 3 distributions with relatively good accuracy (Ref- 
erences 1 and 2) . 

In contrast to the 3 distribution, the glycol-water fluid pumped 
through the porous leading edge is assumed to be distributed at a uniform 
rate over the surface of the porous skin as shown in Figure 3. 

The basic assumption of the proposed prediction method is that the 
minimum glycol flow rate that will still achieve anti-icing (no ice 
accumulated on the leading edge) , results in a freezing temperature of 
the glycol-water mixture, at the point of maximum water catch rate, 
that is equal to the local air temperature. At this point on the wing 
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Figure 3: Distribution of Glycol Flow Rate 

on Porous Leading Edge. 
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(near or at the stagnation point) , the glycol concentration is mini- 
mized; and at the anti-ice threshold the mixture will just begin to 
exhibit the formation of small pieces of ice,. As the glycol flow rate 
is decreased, the extent of ice fermation gradually increases until 
continuous bars are formed spanwise on the leading edge before being 
swept away every few minutes. This latter mode is called natural 
deicing. 

Obviously, to determine the appropriate glycol rate, the local air 
temperature must be known. This temperature will be between stagnation 
temperature and ambient atmospheric temperature. To be conservative, 
one would choose ambient temperature because it will require a lower 
freezing temperature and higher flow rate. However, if a constant flow 
rate is chosen based on the most severe icing conditions anticipated, 
it is probably sufficient to use stagnation temperature as the local 
glycol-water freezing temperature. In this case, the system will per- 
form somewhere between an anti-ice mode and natural deice mode during 
the most severe conditions. Reference 3 shows that the flow rate re- 
quired for natural deicing is only 25% to 50% of that required for 
anti-icing. At le.v.j severe icing conditions, the system will have a 
flow rate equal to or greater than that required for anti-icing. As 
an appropriate compromise, it is suggested that the mean temperature 
between ambient and stagnation be used. 

The method is utilized by following these logical steps. 

1. Find 0 as a function of 
max 

- airfoil shape 

- airspeed 
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- air density 

- droplet diameter 


ORlGINAt PAGIS IS 
OF POOR QUAUT/ 


Any reasonably accurate computer code, such as those discussed 
in References 1 and 2, can be used for this step, 

2, Calculate the water catch rate by the formula 


M = 3 (LWC)V, 
w max * 


( 1 ) 


3. Determine the glycol mass fraction G required to produce a 
solution with a freezing temperature equal to the average be- 
tween ambient static temperature and stagnation temperature, 
using the graph i^hown in Figure 4. For a glycol mass fraction 
between 55 and 80%, no clearly defined freezing temperature 
exists. In practice^ this mass fraction, which corresponds 
with an air temperature below -45 °C, would never be required, 

4. Calculate the fluid flow required to achieve the glycol mass 


fraction G, given a water catch rate M^, by the equation 


GM 


W. = 


w 


f X - G 


( 2 ) 


where X is the initial glycol mass fraction. For most available 
fluids, X is approximately 0.8, 


PREDICTED AMD ACTUAL FLOW RATES 

To test the validity of this method, it was applied to an airfoil 
for which anti-ice fluid rates were determined by tests conducted in the 
NASA Lewis Icing Research Tunnel and reported in Reference 3. 

Those tests utilized an actual airplane wing section. The original 
wing tapered from a NACA 6424215 airfoil at the root (WS 0) to a NACA 
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64j^A4l2 at the tip (WS 216). The wing incorporated increased thickness 
on the forward 30% of the upper surface, a modification proposed by 
R. Hicks of NASA Ames R-sjsearch Center (Reference 4) . The centerline of 
the tunnel v/as at WS 58 of the original wing, where the chord is 63.25 in. 
The predictive method was applied at this station. 

The porous panel mounted on the leading edge consisted of three 
Independently controlled sections, with WS 58 in the center of the center 
section. To obtain the anti-ice flow rates, the upper and lower sections 
were used simultaneously during each run to establish Independent flow 
rate values while the center section was used to obtain minimum flow 
rates for natural deicing. 

The method of determining the glycol flow rate corresponding with 
the anti-ice threshold was as follows. At a given flight condition, the 

rate was set to be well above the anti-ice threshold. The flow rate 
was then reduced in steps, allowing about 30 seconds for the system to 
stabilize at each point, until small flecks of ice began to appear on the 
leading edge in the vicinity of the stagnation point. At the anti- ice 
threshold, the small ice flecks, ranging up to about 3 mm in diameter, 
would form and then be swept downstream in only a few seconds. A glycol 
flow rate lower than the threshold value would allow the ice flecks to 
persist and grow gradually into larger patches before being shed from the 
wing. Complete test details may be found in Reference 3. 

Six test cases were chosen for analysis. Test conditions and flow 
rate data are presented in Table 1. Also shown are the computed values 
of The computer code developed at Ohio State University (Refer- 

ence 2) was utilized. 
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Table 1; Test Data and Predicted Maximum Catch Rate 
for an Anti-ice Porous Leading Edge System 


Case 

''e 

kt 

LWC 

gm/m^ 

d 

ym 

^0 

°F 

a 

deg 

W W 

U 

ml/cm^min ml/cm^min 

p 

max 

I 

96 

1.50 

15 

25 

7.8 

.0255 

.0193 

.208 

II 

96 

2.40 

20 

25 

7.8 

.0436 

.0340 

.314 

III 

96 

1.16 

15 

5 

1.2 

.0420 

.0330 

.219 

IV 

175 

.65 

11 

25 

1.2 

.0240 

.0250 

.259 

V 

175 

.80 

15 

25 

1.2 

.0380 

.0350 

.329 

VI 

175 

1.16 

15 

25 

1.2 

.0540 

.0430 

.329 


To check the validity of the proposed method, the freezing temper- 
ature of the glycol-water mixture at the stagnation point, given P and 
the glycol flow rate, was determined by the following method. 

The water catch rate M was calculated as 


M = p (LWC)V; 
w max * 


then 


XW^ 


G = 


W. + M ’ 
f w 


( 1 ) 


(2a) 


For these tests, X = 0.8. Knowing G, one can determine the freezing 

temperature of the glycol-water mixture using the data shown in Figure 4. 

To determine the sensitivity of the freezing temperature to 

calculations were made for a range of above and below the computed 

value of Pjjjgjjj each case. Also the sensitivity of freezing temperature 

to the glycol flow rate was determined by varying above and below the 

experimentally observed values of at the anti-icing threshold by 20% 

while holding the value of P constant, 

° max 
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HeaulLs of the uflleulutions arc presented in Figure 5 for the. six 
test easea. 

Curves arc presented for the. upper and lower panels, Tha center 
panel, where vras calculated, would be expected to produce, a c,urve 
hetween the upper and lower panel. The consistently lower predicted 
i reeling temperature of the upper panel fluid at the stagnation point 
Is uiMHirently caused by the fact that the calculated d ^ ^ uuderprediets 

nivix 

the actual value for the sharper leading edge of tlie upper panel of 
the tapered wing. The reverse Is true for the lower panel. Tlu're may 
also he variations In the I.WC across the model spun. 

In every ease o.xco,pt for Cose 1, the midpoint between the upper 
and lower panel curves predicts a free/.lng temperature at the point 
of that lies between the ambient and total air temperature, a 
result to be expected. This caleulated freezing temperature is based 
on the predicted value of bbe observed glycol flow rate, the wind 

tunnel test conditions, and tlie properties of a glycol-water solution. 

The sensitivity of the predicted freezing temperature to errors 
in the value of f and to variatiovts in the glycol flow rate are illus- 
trated for each case. 

It 1? "ow useful to compare, actual values of Wj^ with those pre- 
dieted by the method presented heroin. Ivcsults are shown in Table 2. 
Note that the glycol flow rate is converted from volume flow rate to 
mass How rate by the specific gravity of the original fluid, which 
is 1 . I . 
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kt 

CASE ( Tq-25 

LWC = iS gm/m^ 
<< - 7,8 dog. 




Figure 5: Effect of on Local Fluid Freezing Temperature and 

Effect of Flow Rate on Local Freezing Temperature at 

Predicted S • Case t. 
max 
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LWC - 2.4gm/m^ 
^=7.8 d«g. 




Figure 5: Effect of Local Fluid Freezing Temperature, and 

Effect of Flow Rate on Local Freezing Temerature at 

Predicted 3 . Case II. 

max 
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CASE /// 


Vg = 96 kt 
To= S -F 
LWC = U6 gm/m^ 
- L 20 deg. 

d “ !5/<m 


•0 “ predicted 3 
for test Max 

conditions 


Test Wf” 
values used^ 


0 . 15 / 0-20 


Calculated ^ Max ^^Max'"'^^^ 
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Figure 5; Effect of 3 on Local Fluid Freezing Temperature, and 
° max 

Effect of Flow Rate on Local Freezing Temperature at 


Predicted 3 


Case III. 
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Figure. 5: 


FT feet of on Local Fluid Firoealng Temperature,, .and 

Ef,Pect of Flow Rato on Local Freosing Temporature at 
Frotlicted Case V. 
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Table 2: Comparison of Predicted with Actual Test Flow Rates 

at Anti-icing Threshold 


Case 

3 

max 

gm/cm^-min 

T + T 
a 0 

predicted 

gm/cm^min 

W- + W 
U L 

Prediction 

2 

deg F 

2 

gra/ cm^-min 

Accuracy 

% 

I 

.208 

.0837 

23.9 

.0179 

.0246 

-27.2 

11 

.314 

.202 

23.9 

.0432 

.0427 

+ 1.2 

III 

.219 

.0666 

3.9 

.0415 

.0413 

+ .5 

IV 

.259 

.0878 

21.5 

.0228 

.0269 

-15.2 

V 

.329 

.137 

21.5 

.0356 

.0401 

-11.2 

VI 

.329 

.199 

21.5 

.0517 

.0533 

- 3.0 


These results show that the method of prediction of anti-ice 
threshold flow rates presented herein predicts flow rates with an 
average error of less than 10 percent of the experimentally determined 
flow rates. This is believed to be excellent, considering the estimated 
accuracy of the prediction and the fact that the anti-icing 

threshold characteristics tended to persist over a relatively wide range 
of values of glycol flow rate, making it difficult to obtain fine reso- 
lution of the anti-icing threshold. Furthermore, the determination of 
the anti-icing threshold was somewhat subjective, depending on the 
judgement of the experimenter. 

Therefore, this method appears to provide a reasonably accurate 
determination of the flow rate required to assure continuous anti-icing 
performance at a given icing flight condition. For less severe icing, 
an excess of anti-ice fluid is available at the leading edge. For 
more severe icing, there will be a gradual transition to a natural 
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deicing mode which generally provides adequate protection aga.inst ice 
buildup in flight, although this would have to be evaluated for each 
particular alrc&'aft. 

It is recommended that additional testing be done on different 
airfoils under a variety of test conditions to verify the generality 
of the method presented in this report. 
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